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Abstract 

We study counting statistics of electric current pumped by pulses of an exter- 
nal field. The fluctuations depend on the pulse shape, and can be minimized 
by choosing the pulse shape properly. For an optimal pulse shape, the fluctu- 
ations are reduced to the dc level, i.e., they do not depend on the duty cycle 
of the signal. We develop an approach that allows to calculate all counting 
statistics for various driving fields, optimal and non-optimal. The statistics 
depend in an interesting way on the analytic structure of the field time de- 
pendence, and display an analogy with coherent states and instantons. 

PACS numbers: 72.10.Bg, 73.50.Fq, 73.50.Td 
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I. INTRODUCTION 



One of the fundamental problems of "single electronics" is that of quantized charge 
transfer monitored by an external field [Q . A real device such as electronic pump or turnstile 
is operated by a periodically alternating field which drives electric current at a rate of one 
electron per cycle 0. Due to a microscopic scale of the currently studied systems, the issue 
of current fluctuations, classical and quantum, becomes very important. Ideally, in order to 
reduce current fluctuations well below the level of the ac signal, the pump should transfer 
charge adiabatically For a perfectly operating pump, the cycle duration T is much longer 
than h/Ec, where Ec is charging energy. If the adiabatic limit E^T ^ ^ is achieved, the 
current fluctuations are exponentially small. 

In practice, the adiabatic parameter can be of the order of one which makes the analysis 
very difficult, since one has to consider many effects simultaneously: co-tunneling 0, finite 
relaxation rates quantum fiuctuation of charge ||^, etc. In this paper we consider the 
problem in the non-adiabatic limit E^ <^ fi/T, where the charging energy can be ignored. 
With the usual assumption that Fermi-liquid quasiparticles' interaction vanishes near Fermi 
surface, the pump can be described by single particle scattering amplitudes periodically 
varying in time. In this case the only effective mechanism of the current fluctuation sup- 
pression is due to Fermi statistics which makes electron transmission events correlated in the 
time domain. The effect of Fermi statistics on the dc current fluctuations has been studied 
0. We shall extend this theory to the ac current. 

From the point of view of the pump quality, the extreme non-adiabatic limit we consider 
is the least efficiency mode. So, it gives an upper estimate of the fluctuations at low temper- 
ature, and provides a reference for understanding real systems where non-adiabatic effects 
are strong. Actually, for the non-interacting fermions one can develop a complete theory 
which gives not just the mean square of current fluctuations, but all statistics of transmitted 
charge [|l^ . Although the free fermion counting statistics problem captures only part of the 
physics relevant for operation of a real device, it is interesting enough and non-trivial by 
itself. 

It turns out that the fluctuations strongly depend on the pulse shape of the driving ac 
signal. The character of this dependence resembles H-theorem. It was shown recently by 
two of us that at flxed average current the fluctuations level is bounded from below by 
the fluctuations of the dc current with the same mean value Moreover, there exist 

"optimal" ac signals for which the minimum is reached, and the fluctuations remain on the 
dc level (see Fig. 1). It is interesting that for optimal signals the fluctuations are independent 
on the relative pulse width r/T (the signal duty cycle), no matter how sharp the pulses are. 

The case is appealing to an analogy with coherent states that minimize quantum me- 
chanical uncertainty. Similarly, the optimal ac signals drive current in such a way that 
quantum noise is reduced to the dc minimum level. Coherent states are known to possess an 



interesting analytic structure [T^, and one encounters a similar situation in this problem. 
We study analytic structure of the problem, and show that it is related to the group of 
modular transformations. The analytic structure facilitates the study of counting statistics. 
We develop a method discussed in our previous work and use it to calculate complete 
counting statistics for various interesting examples of the driving signal. 

To be speciflc, we study a one dimensional model of a Fermi gas of electrons transferred 



2 



through a potential barrier by an external ac field. Electrons are incident from the left 
and the right reservoirs which supply zero temperature Fermi distribution. In the ac field, 
the scattering becomes non-elastic, and electrons with different energies interfere due to 
their Fermi-statistics. We treat the problem as a multichannel scattering problem with the 
scattering amplitudes given by the Fourier components of the external field, and present a 
general formalism that gives counting statistics. 

We compute the distribution exactly for a particular class of external fields, which in- 
cludes the optimal signals. For these fields the law of time dependence is a periodic analytic 
function of time given by a rational functions of the "circular" variable z = e*^*. We find 
exact statistics of the charge transfer, and show that it displayes interesting features ex- 
plained by the analytic character of the time-dependence law. In particular, we shall discuss 
Lorentzian pulses of voltage with quantized fiux: c / V(t)dt = n$o, where n is an integer, 
and $0 = " is the fiux quantum. We find that the pulses represent a quantum analog of 
the classical picture of n independent attempts to transmit electrons through the scatterer. 
We arrive at this result by means of the method of [T^ combined with a special treat- 
ment which allows one to reduce the problem in an infinite- dimensional space to a finite 
dimensional problem (cf. Ref. [0). 

The paper is organized as follows. In sections 2 and 3 we describe the model and state the 
mathematical problem to solve. We introduce a generating function of the charge transfer 
statistics and express it as the determinant of a scattering operator in an infinite-dimensional 
space. For computing the determinant it is important to find physically meaningful regular- 
ization. The regularization problem is treated in section 4. We find that the regularization is 
sensitive to gauge transformation of the electromagnetic field, and choose the gauge so that 
the regularization becomes simple. In section 5 we address the question of noise minimiza- 
tion and recall that the signals that give minimal noise have particular analytic structure: 
rational functions of the circular variable e'^* analytic inside the unit disc 0. We show 
how this result follows from our expression for the generating function. Next, in section 6 
we treat the case of a rational time-dependence of the external field for which we reduce 
the problem to finding the determinant of a finite-dimensional matrix, and compute the 
probability distribution explicitly. Section 7 contains the discussion of the symmetries of 
the problem. We observe that the system possesses a symmetry group of conformal transfor- 
mations isomorphic to PSL{2,IV). Section 8 treats several interesting examples of driving 
signal which reveal interesting features of the counting statistics. In section 7 we summarize 
our discussion. In order to make the readers' burden less heavy, some technical details are 
moved from the main body of the paper to the Appendices A and D. For completeness, we 
review relevant results of 0] and \T2] in the Appendices B and C. In Appendix E we recall 
gauge transformations that relate the ac fiux problem to the ac bias problem, which are 
essential for the issue of the determinant regularization. 



II. MODEL 

We consider the following model of a microscopic contact. One-dimensional ideal Fermi 
5as scatters off a potential barrier U{x), so that the Hamiltonian of the system is 
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H = ^^{c 



.d_ 

dx 



-a{x,t)\ + ev{x,t) + U{x) 



(1) 



where and are the canonical operators of electrons, U{x) is the scattering po- 

tential, a{x,t) and v{x,t) are the external electromagnetic field vector and scalar potentials, 
respectively. We shall treat electrons as spinless and non-interacting. The scattering poten- 
tial U{x) will be described by the matrix of scattering amplitudes 



A 



All 
Alr 



Arl 
Arr 



(2) 



so that the scattering states have the standard asymptotic form: 




Allc-'"'', 
Alrc''^^ 

Arlc-'''^ 
Arrc'^^ 



X - 

X - 

X 

X - 



— oo; 
+00; 

► —00; 
+00. 



(3) 



We would like to study the response of the system to a time-dependent external field, 
electric or magnetic. Let us recall that by a gauge transformation one can go from a problem 
with the electric potential v{x, t) to that with the vector potential a{x, t) and vice versa (see 
Appendix E). In the limit of instant scattering which we shall assume from now on, the 
external field may be taken into account by introducing a time-dependent phase in the 
scattering amplitudes: 



A ^ k{t) 



All 

Alrc^'^^'^ 



ArlC 
Arr 



(4) 



"Instant scattering" means that the external fields a{x,t) and v{x,t) (and thus e**^*^*^ as 
well) vary slowly in comparison with the scattering time Tgc ~ h\dAa/3/dE\, where Aa/3 are 
the scattering amplitudes, E is the energy of electrons. The physical meaning of r^c is the 
time that the particle spends inside the scatterer. Technically, neglecting Tsc means that the 
scattering amplitudes Aap do not depend on the energy. 



III. STATING THE PROBLEM 

Now, after the model is defined, we are ready to state the problem of the counting 
statistics of electron transmission. Namely, we are interested in finding the probabilities of 
a given charge transfer over a fixed time interval. We assume that the counting time T*^°) is 
much longer than the period T of the external field. Let us observe the system under the 
action of the external field during a long time interval T'^°^ The field from now on is treated 
only as the phase factor e**^*^*^ in Eq.(^). Let us denote the probabilities of transmitting 
exactly n electrons during this time by -P^°^. The probabilities can be conveniently combined 
into a generating function 

+00 

X^''W= E Pi'^e^'""- (5) 
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This function of the auxihary parameter A will be the main quantity we shall work with. It 
encodes all the information about the statistics of the charge transfer. The moments of the 
distribution are given by the coefficients of the Taylor expansion of x^'^K^) around A = 0. 
In particular, the average charge transfer is 



W = -^1; (6) 



A=0' 



and the dispersion 



{{n')) = {n') - in)' = '^^L,^^ X^'\>^) ■ (7) 



Notice that even if we apply no external field there exist equilibrium fluctuations. At zero 
temperature {{n'^)) is growing with time as InT'^'') |]l3i. The zero temperature equilibrium 



fluctuations can be neglected if the external fleld e*'^'-*^ is periodic in time: in this case the 
fleld-induced fluctuations grow linearly with T^'^^ and dominate over the equilibrium noise 



IJ] . In this paper we disregard the equilibrium noise and study only the effect of the periodic 
fleld, and only briefly consider the non-periodic case and the equilibrium noise in Appendix 
D. 

There exists a useful description of the field-induced statistics in terms of the quantities 
that do not depend on the counting time. Let us introduce 

Pn = {probability of transmitting n electrons per cycle}. (8) 

The advantage of P„ is that the total generating function x^^H-^) ^oi a long time interval 
2"(o) ^ x can be written as an exponent: 

x(°)(A)~[x(A)r'/^ as T(°)^oo, (9) 

where 

+00 

X(A)= Y: Pne''"- (10) 

71= — oo 

is the generating function per one cycle of the field. 

This factorization can be explained by recalling that the zero temperature equilibrium 
fluctuations are essentially due to the counting beginning and ending. It is completely anal- 
ogous to the fluctuations in one-dimensional ideal Fermi gas (here the role of the dimension 
is taken by the time). One can regard the fleld-induced fluctuations as "extensive in the 
time domain", and the equilibrium fluctuations as "boundary effects". As usual, to treat 
the "bulk" effects separately from "surface" ones, we adopt periodic boundary conditions, 
which leads to Eq.(^. 

More formally, let us flrst suppose that the phase factor is periodic: e*'^*^*"''"^^ = e*'^*-*-'. 
Then the factorization (^ can be understood by the following argument. We observe the 
system during a large number of fleld cycles and impose boundary conditions periodic in 
time. The spectrum of energies then is discrete with the spacing hQ/N = h/NT. Scattering 
is possible only between energy levels separated by a multiple of HQ = h/T. Thus, we have 
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non-interfering copies of the scattering problem with discrete energy levels spaced by HQ. 
omce N = T^^)/T, this proves our formula (y). 

If the phase factor e*'''^*-' is not periodic, the factorization still holds. In this case, since 
the voltage ^0 is periodic, the phase factor is quasi-periodic. The generating function of 
the statistics per cycle is constructed in the following way. One has to find the smallest 
positive e for which the phase factor e*'^(*)+*''* is periodic. Then x(A) = Xi"^('^)X2(-^)5 where 
Xi and X2 are found for the periodic phase factors e*v(*)+*'^* and e«¥'(*)+«(<:-2'r)t^ Pq2~ ^]-^g proof 



and discussion of the quasi-periodic case we refer to our previous paper |T^. In this paper 
only periodic phase factors will appear. 

The above argument shows that the field-driven contribution x{^) can be extracted 
by closing the time axis into a circle of period T, thus quantizing the energy with the 
quantum hVt = h/T. The problem becomes a multi-channel scattering problem, where the 
channels represent the discrete energy levels. For such a problem xW can be quite generally 



expressed in terms of scattering amplitudes [12]. In this method the function xW depends 



on a vector argument A = (Ai, . . . , Aat), being the number of channels. Fourier transform 
of xW gives the probabilities of scattering between channels. If A is the matrix of the 
multi-channel scattering amplitudes, n is the occupation number operator (diagonal in the 
energy representation), A = Diag{e^'^'^ , . . . , e*'*'^), then 

x(A)=det(l + n(S-l)), (11) 

where 

S = AA. (12) 

The proof of this formula is reviewed in Appendix A. 

In our treatment of the periodic problem the operators n. A, A, etc. act in the linear 
space 

H = V®V (13) 

of left and right states with discrete energies. Formally, V is the C°° space of the states of the 
discrete spectrum, and by Fourier transformation it can be treated as the space of periodic 
functions of time, with the period T. Two copies of V correspond to the left and right 
channels. Further the computations will be performed in the basis (a, m), where a G {-R, L} 
specifies the side, m G Z labels the Fourier harmonics. To be specific about notation, we 
shall sometimes write operators F acting in H in the form 



Flr Frr 

where Fap are operators acting in = C°°. 

In the application to our problem, the scattering matrix is time dependent: 

l_( B A*[f{z)Y 
\Af{z) -B* 



(14) 



(15) 



where A and B are the transmission and refiection amplitudes for the potential barrier 
U{x), -|- = 1] z = exp(27rzt/T) is the variable on the unit circle; the phase factor 
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f{z) = exp{iip{t)). We choose the formal variables Aj to count only total charge transfer, 
regardless of the energy of transmitted electrons: 



^L,m — 0, ^R,m — A . (16) 

Then from Eq.(|l2D we get 

where the functions f{z) and [/(-z)]* = f{z) should be understood as operators acting in 
V by multiplication: 

^m\f{z)\n)= i^J^. (18) 
We consider the system at zero temperature, thus 

n a,m) = i ' ' (19) 

' ' 10, m > 0. ^ ^ 

By that, formally, all operators in ( [Tl| ) and (|I2|) are specified, and we can proceed with the 
calculation. 



IV. REGULARIZATION OF THE DETERMINANT 

Now, the problem is to find the determinant (pJ] ) of an infinite matrix. Because of the 
infinite dimensionality, this determinant needs to be understood properly. Eventually, to 
compute the determinant, we are going just to cut a finite submatrix of 1 + n(S — 1) at some 
positive and negative energies way off the Fermi level, and to take its determinant. However, 
one has to be careful and to make sure that the infinite parts of the matrix thrown away 
do not matter. Physically, the reason is that all incident states with very high energies are 
empty, and remain such after scattering off an alternating field. Similarly, all states with 
very low energies are deep in the Fermi sea, and remain doubly occupied throughout the 
scattering. Thus, in order that our regularization procedure has no effect on the determinant, 
we have to assure that the unitarity of scattering is preserved when a submatrix is cut. Then 
the infinite submatrices that we remove will contain only the contribution of "Fermi vacuum" 
states, and will not affect the counting statistics, and thus the determinant. 

More formally, we have an infinite number of states below the Fermi level in the left and 
the right reservoirs to be put in a one to one correspondence. Shifting states in one reservoir 
with respect to the states in the other one can be achieved by a gauge transformations 
different for the left and the right reservoirs. The key observation is that the matrix 1 + 
n(S — 1) is not invariant under a gauge transformation. (The transformation rule of the 
particle number operator n is reviewed in Appendix E, and the matrices A and S transform 
the same way as the density matrix.) Of course, since the determinant of 1 + n(S — 1) gives 
counting statistics, the determinant regularization must depend on the gauge transformation 
in such a way that the regularized determinant is gauge invariant. 
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In order to clarify the relation between the regularization and the gauge transformations, 
let us consider a system with no barrier. Then A = 1, B = 0, and the scattering is 
only forward, no backward. Then the matrix S given by (|1^) becomes diagonal and time- 
independent: 



e-'^ 




AX I 5 



(20) 



In this case, one could try to compute the determinant of 1 + n(S — 1) by using "naive" 
regularization, i.e., by simply cutting all columns and rows above some large positive and 
below some large negative energy. Then xW = det(S) = 1: no transport for any v^(t), 
which is a meaningless result. The problem becomes even more striking if one thinks of a 
gauge transformation. Under a gauge transformation the S given by (|20|) does not change, 
and n is transformed according to the rule (^) (see Appendix E). For example, if the gauge 
phase 0(t) = nQt, in the energy representation one gets 

n'^{E) = ULiE) , n'^(E) = n^(E - nhQ) . (21) 

Now, by using the naive regularization one finds 

X(A) = e"^ , (22) 

which means that the result is not gauge invariant. 

At this point one can conclude that the correct regularization must change together with 
the gauge. For the general scattering problem (|15D we proceed in the following way. We 
choose a gauge transformation so that the scattering matrix A becomes time independent. 
Then S is also time independent, and all phase factors are shifted to n. We argue that 
after such gauge transformation one can use the naive regularization. The transformation is 
chosen so that all dependence on the ac field is moved to the occupation of the incoming one-, 
two-, or many-particle states. The advantage of going to a purely elastic scattering is that 
all outgoing channels that enter the scattering unitarity relation will have equal energies. 
Thus, while removing from 1 -|- n(S — 1) all states with energies below some large negative 
energy we preserve the unitarity of scattering. In other words, the states that interfere at the 
scattering are either both included in the truncated matrix, or both are removed. Therefore, 
for our gauge transformation the naive regularization is meaningful. Moreover, it is clear 



that such transformation is unique, unless S is diagonal |jT5[ . 

On the basis of this consideration, we argue that the determinant (|1^) must be under- 
stood in the following way: 

X(A) = "det"(l + n(S-l)) = "det"(l + n(^[^^^)]* ° ) (S - 1) ^[^"^^ ° 

= det(l+(^|;^) J)n(f^^Q)]* °)(S-l))=det(l + n(S-l)), (23) 



where 



|^|2(^g-iA _ 1) + 1 A*B*{e'^ - 1) 
-AB(e-'^ - 1) \A\^(e'^ - 1) + 1 



(24) 
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is the time-independent scattering matrix, and 



V i; V 1 



(25) 



is the time-dependent occupation number operator. In the last hne of (^) we can assign 
an unambiguous meaning to the determinant (cf. Jill). Namely, we note that S is unitary 
and has unit determinant. Also, n tends to at high energies and to 1 at low energies. 
Therefore, the matrix (1 + n(S — 1)) behaves at infinity like a block-diagonal matrix, with 
2x2 blocks each having unit determinant. We define the determinant in Eq.(PBD by cutting 
the matrix along one of these blocks at infinity. In such a way the determinant is well defined 
and depends essentially on the matrix entries around the Fermi level. Below we shall be 
able to compute it explicitly for a special choice of f{z). 



V. NOISE MINIMIZATION 

It turns out that the analytic structure of f{z) plays an important role in the current 
fluctuations. In particular, the functions f{z) which minimize the noise {(n^)) for a given 
average charge transfer (n) belong to the class of rational functions. The problem of noise 
minimization has been considered in detail in |P, and here we briefly review the result. The 
rate of charge transfer {n) is simply proportional to the phase gain per period: 

{n) = \A\'^, A^ = ^Tgf{z)\lzl (26) 
while the noise {{n^)) depends on the whole function f{z): 

For completeness, both expressions are derived in Appendix B. 

The variational problem arises of finding the function f{z) which defines a map of a fixed 
degree N = A(y9/27r of the unit circle into itself and minimizes the noise functional (0). In 
Appendix C we review the proof of Ref. IQ that optimal f{z) is analytic either inside or 
outside the unit circle \z\ = 1. In other words, its Laurent expansion f{z) = J2n=-oo '^n^"' 
contains either only non-positive or only non-negative powers. Such a function can be 
written as 

/w = ni^. (28) 

i=i "i^ 

where either all |aj| > 1, or all |aj| < 1. The corresponding time dependence of the phase is: 

/ X 1 / (1 - |a^^sin^](^-^,) \ , 

ip{t) = 5:tan-i \ , ' ^' \, , + 00 , (29) 



i=l 



cos^l{t — tj) — 2\ai\ 



where = argaj/f2, 0o = Z^jargc^i- Thus the optimal phase time dependence is a sum of 
N "elementary excitations", or "kinks", each corresponding to a 2tt phase change of the 
scattering amplitude per one cycle of the signal. 
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For any of such functions the mean square fluctuation of the transmitted charge (((en)^)) 
is equal to e'^\AB\'^N per cycle. It is remarkable that the noise does not depend on relative 
displacement of the kinks in the time domain, as well as on their durations. The degeneracy 
is described by 2N — 1 real parameters. 

It is interesting that for the time dependence ( pH]) all the probabilities P„ can be computed 
and admit a simple interpretation of non-interfering attempts of electrons to pass through 
the barrier (see Examples 1 and 2 below). In fact, the class of functions to which our method 
applies is broader: it includes all rational functions (^) regardless of the location of the 
poles. (The phase time dependence then has the form ( PP] ) with arbitrary signs of different 
terms.) We show that for any such f{z) the generating function xW can be expressed 
as the determinant of a finite matrix, and that it yields only a finite number of non-zero 
probabilities Pn- 



VI. COMPUTATION 

Let f{z) be a rational function that maps the unit circle into itself: \f{z)\ = 1 for \z\ = 1. 



Then it has the form (p8|) with arbitrary \ai\ ^ 1, not necessarily all inside or all outside the 
unit circle. Let 

P{z) = Uli{l-a;z), (30) 
so that f{z) = Q{z)/ P{z). We shall also use the functions 



1 



kez 



p-\^) = j^, = j:cr^\ (31) 

where the Laurent expansions are chosen to converge on the unit circle. We shall treat the 
functions (|30D and (|3l| ) as operators acting in V by multiplication. 
Now we use the following trick to compute the determinant. 



X 



(A)=det(l + n(S-l)) =det(l+(^/^ ?)^(^0^ J)(S-1)) 



-det(i+(^; g-OKo g)(^-^))- (32) 

Here we performed the gauge transformation in both left and right channels simultaneously. 
This does not change the determinant. 

Simple computations show that the matrices P~^nP and Q~^nQ have the following form: 
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p-^nP 



Q-^nQ 




(33) 



(34) 



blank spaces stand for zeroes, crosses mark the Fermi level, asterisks denote arbitrary entries 
not used in calculations. Therefore 



X(A) = det 



5n 



* 


* 


* 


* 


* 


* 








Xni 




Xnn 


* 


* 


* 


* 


* 


* 



det 




XiN 

Xnn 



where 



Xij 



(S-J) 











+ 16, 



(35) 



(36) 



are 2x2 matrices. 

The determinant (|35|) is finite and determines xW ^ function of the parameters {aj}, 
A, and B. As a function of e^^ and e~^^, the determinant is a finite degree polynomial, since 
all entries of Xij are such (see Eg. (|3^) ) . As a result, there is only a finite number of non-zero 
probabilities Pk in the expansion ([To| ) for x('^)- 



VII. THE PSL{2, R) SYMMETRY 

Before we turn to the discussion of examples, let us mention that the system possesses 
an interesting symmetry group PSL{2, R) defined as the group of real unimodular matrices. 
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The symmetries are realized by the group G of hnear-fractional transformations preserv- 
ing the unit disc 1^1 < 1 (this group is 3-dimensional and isomorphic to PSL{2,'R)). We 
claim that: 

If two phase factors f{z) and f{z) are related by such a transformation {f{z) = 
f{g{z)) for some g G G), then the generating function of the distribution x(A) is 
the same for f{z) and f{z). 
Indeed, let us notice that xW is expressed as the determinant ( pUD of an operator 
in the space of functions on the circle \z\ = 1. Also, any g & G commutes with the 
occupation number operator n at zero temperature [g does not mix positive and negative 
Fourier harmonics). If we perform the conjugation by g, the expression ([TT| ) remains the 
same, except that f{z) gets replaced by f{z). The determinant however does not change 
under a conjugation. This proves our statement. 

Let us note that f{z) is uniquely determined by its zeroes (or, equivalently, by its poles 
1/a*). The sets of the parameters for f{z) and f{z) are also related by g: ai = g{di). 
Therefore, as a function of {ai}, xW is invariant under the simultaneous mapping of all 
hj geG. 



VIII. EXAMPLES 

Let us illustrate our discussion by actual computing the probability distribution according 
to Eq.(|35D for several specially chosen functions (p^). 

Example 1. The simplest case is = 1. Obviously the transformation h-* l/oj just 
switches the direction of the charge transfer. Therefore, without loss of generality let \a\ > 1. 
Then by using the expansions (|3lD one gets 

9ii = 1, 
Pn = 0, 

X(A) = |A|V^+ (37) 

Thus we have exactly one attempt to pass the barrier per period with the probabilities \A\'^ 
to pass and to rebound. Such a situation already occurred in the problem with constant 
voltage |I12|. Constant voltage is a special case of this example corresponding to a = or 
a = CO (i.e., / = exp{±iQt) ). 

Example 2. Let now A^ > 1, and the Oj, i = 1, A^, be all inside or all outside the unit 
disc (see Fig. 2. a). Again, without loss of generality let all \ai\ > 1. Now, we obtain 

Pij = 0, 

X(A) = (|A|V^ + l^n^. (38) 

We see that xW contains A^ equal factors, each corresponding to one of the factors of f{z). 
This means that each "elementary excitation" in ( ^9]) corresponds to one attempt to pass the 
barrier with the one-particle outcome probabilities. Of course, the actual scattering state in 
this well as in Example 1, is a many particle coherent state. One could expect this 

result when all Oj = oo (this again corresponds to constant voltage). However, here we get 
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the same expression for the superposition of any set of "elementary excitations" (Example 
1) of the same polarization with no dependence on the values of Oj. It is a surprising result. 
Of course, in this example different sets of {oj} cannot be transformed into each other by the 
symmetry group G, therefore we cannot explain this invariance merely by the PSL{2, R) 
symmetry discussed previously. At present we look at this invariance as at a miracle and 
admit this to be a consequence of a broader group of symmetries. 

We shall also list the answers for two other examples of f{z) which demonstrate the 
interference of "excitations" with opposite polarizations. 

Example 3. N = 2. \ai\ > 1, \a2\ < 1. Then one has 



al — al \a1a2 —02 
1 / ai —a\a2 



where 



X(A) = 1 -2F + F(e*'' + e-*^), (39) 

F=iAni? H^~"^"f . (40) 

10-1—0,21 

Note that in this example the probabilities of the electron transfer in both directions 
are equal, and thus average charge flux is zero, although for arbitrary parameters ai and 
02 the signal can be asymmetric (see Fig. 2.b). In this case x(-^) cannot be represented 
as a product of independent contributions of the two "elementary excitations", but exhibit 
their interference. To make this clear, let us show what happens if one tries to factor the 
generating function: 

;^(A) = {u + we'^){u + we-'^) , (41) 

where the "probabilities" u,w = |(1 ± — 4-F)), so in this case there is no natural relation 
between the factors of xW and of f{z). 

Example 4- N > 1. ai = a, = b for i > 1. \a\ > 1, \b\ < 1. Computations can be most 
easily performed for the case 6 = (which corresponds to a constant voltage applied against 
one "elementary excitation"). Then by using the PSL{2,TV) symmetry one can extend the 
result to arbitrary b: 

X(A) = (1 - 2F + F{e'^ + e-'^)){\A\^e'^ + \B\Y~\ (42) 

where 

^=i^n^i'n— (43) 

\a — op 

Like in Example 3, there is an interference of the kinks of opposite sign. It is interesting 
that the factors of x(A) can be interpreted by saying that one positive and one negative kink 
interfere and form a "neutral" system with which other kinks do not interfere. The degree 
of interference is measured by A = |1 — a*6|/|a — 6|, and varies from to 1 depending on 
how much the kinks overlap in time, or by how much their durations differ. For example, 
A ^ 1 if the kinks almost not overlap, and then x(A) factors into separate contributions 
of independent kinks: the "probabilities" u and v in (|41|) become just the one particle 
probabilities \A\'^ and 
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IX. DISCUSSION 



To summarize, we found the probability distribution for the charge transfer under the 
action of a periodic external field described by a rational function of z = e*^*. According to 
our previous remarks about the equivalence between magnetic and electric fields (also, see 
Appendix E), we can treat the field f{z) of the form (pH]) as the alternating voltage V{t): 



1 - 




2 


1 + 


\ai\ 


2 - {a*z + aiZ-^) 



p1 — 2pi cos(f2t — (fi 



(44) 



where = pie^'^\ 

This expression represents V{t) as a sum of "elementary excitations". Each elementary 
excitation alone represents one attempt of an electron to pass the barrier, with the prob- 
ability to pass given by the transmission coefficient of the barrier. We found that if the 
excitations are of the same sign, then these attempts do not interfere. The general formula 
(^) holds for any V{t) composed of elementary excitations of arbitrary signs. It predicts a 
nontrivial interference of the excitations of opposite signs. 

Each elementary excitation in Eq. (H3) can be written as a sum of Lorentzian pulses: 



with the width = ln(p^^), centered at tk = fi^Vfc- Each of the pulses carries flux 
proportional to its area. The flux is quantized: 

Q 2T]fdt he ^ 

2 = - = "^0 . (46) 



e {t - tk - rriTy + 

The flux quantization is just another way to say that each elementary excitation corresponds 
to a 2tt phase shift. 

Using the representation ( ^5] ) our results can be to some extent translated to the non- 
periodic case. The limiting form of the excitation as the period T oo {Q 0) is 

a single Lorentzian pulse with the area JV{t)dt = h/e. If many such pulses all of the 
same sign are generated continuously at a finite rate so that average current dominates 
over the equilibrium fiuctuations, then the calculation will tell that each pulse corresponds 
to a "one-electron like" attempt to pass the barrier, and the distribution of outcomes is 
binomial (exactly as for constant voltage). The interference of the field-driven current with 
the equilibrium noise will be discussed elsewhere. 
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X. CONCLUSION 



We studied quantum counting statistics of an ac current driven by pulses of external 
field. There are special pulse configurations which create many particle coherent scattering 
states in which the quantum noise is reduced to a dc minimum. The analytic structure 
of such states is studied, and put in connection with the modular symmetry group of the 
problem. A general method to calculate counting statistics is presented and applied to the 
coherent states and to other states "naturally related" to them. The counting statistics are 
found to be binomial and "generalized binomial" , respectively. 
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APPENDIX A: EXPRESSION FOR THE GENERATING FUNCTION 



In this Appendix we review the proof of (2.f) from |jT2[. First we express the determinant 
(2.f) as 

x(A)=det(l-n + nS)= S^Z IlO- ' ^i) II (Al) 

{ii,-,ik} «=«a 

where the summation is over all subsets of channels {ii, . . . ,ik}. Here Sll[['J^ denotes the 
determinant of the submatrix of S formed by the entries in the rows {zi, . . . , z^.} and in the 
columns {ji, . . . , j^}. From (|12D it follows that 

StZ= E e^(^-+-+^^.-^H--A.J|A^-^|2^ (A2) 
{iivJfc} 

where the determinants 

^n::t = E^(^X''' (A3) 
p 

are n— particle scattering amplitudes. Note that 

are the probabilities of the many-electron scattering from the channels ii,...,ik to the 
channels ji , . . . , . This proves that 



X(A) = E ^n.....b......e^^^-+-+'--'---'-^ (A5) 

{il,...,if,} 
{il.---.Jfe} 

is the generating function for the probability distribution Pii...it\ji...jk charge transfer. 
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APPENDIX B: EXPRESSIONS FOR {N) AND ((iV^ 



From (^) we can derive expressions for the total charge transfer and its dispersion in 
terms of the external field e*'^*-*^ We shall perform the calculations for the periodic field 



/(e*™). One can write 
d 



{n) = —i 



d 



d 



A=o ' dX 



A=0 



;S - 1)) = Tr 



n 



-AB 



A=0 

A*B* 



:i + n(S-l)) 



In the time representation [z = e 



int\ 



n<0 



n Sn 



1 



Z2 



where 6 is infinitesimal positive. Thus one gets 



Z2 



Z2{l + S)~Zl 





Z2{l + 5)-Zi 





1 



z2{l+5)-zi 



Therefore, 



(n) = \A\ 



lim 

Z2-*Zi 



f{z,)nZ2)-l\ dZi 



Zl - Z2 



2m 



where A(f is the phase change per period. 
Similarly, 



-Tr 



A=0 



^ d 



92 

Inx(A) = -Tr^ 
\ 2 



A=0 



ln(l + (fi(S - 1)) 



A=0' 



fs-r 



n- 



Tr 



-AB 



n 

dzi 

lim 

27rz 22^2^1 



M|2 



\A\ 



n 



-\A\^ 
-AB 



A=0' 



-A*B* 
UP 



n 



\AB 



dZ3 

2m 



{\A\ 



f{z^)r{z2) + i 



+ \A\'\B\ 



2\Ar\B\ 



{Zi-Zs{l+6)){ZS-Z2{l + S)) 

2UP 



-A*B* 
UP 



f{z,)nzs) + f{z,)nz2) 

{zi - zsil + 6)){z3 - Z2{1 + S))^ Zl - 2:2(1 + 6) 
dzidz2 f{zi)f*{z2) - 1 



{2m) 



[Zl - Z2) 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



The last line shows that {{n'^)) depends in a non-trivial way on the whole function f{z), 
unlike (n) which depends only on the total phase shift Aip per period. This recovers the 
result of for a periodically varying field. 
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APPENDIX C: VARIATIONAL PROBLEM 



In this Appendix we review the proof of ^ that the variational problem of minimizing 
((n^)) for a fixed value of (n) is equivalent to the analyticity of f{z) either inside or outside 
the unit circle. 

We decompose f{z) into a sum of f+(z) and f-{z) which are analytic inside and outside 
the unit circle, respectively, 



oo 



n=0 

oo 



n=0 



Then, by Cauchy theorem. 



i|A(.W!W = /^/_W9/;W = o. (C2) 



Therefore, 



and 



{n) = -\A\'j^J{z)dr{z) = -\A\'j^[U{z)dfl{z) + Uz)dnz)] 

= \A?T.<\<\"-\<^nW (C3) 



dz 



((n^)) = 2\A\'^\B\ r r , 

\A\^Wj^^(dU{z)-df^{z))nz) 
= \A\'\B\' j ^^[dU{z)fl{z) - dUz)nz)) 

= \A?\B\'Y.<\<? + K\% m 

n 

By comparing the two expressions we see that at fixed {n) the fluctuation ((ri^)) is minimal 
when /+ or /_ vanishes: then {{n'^))min = \B\'^\{n)\. 



APPENDIX D: NON-PERIODIC SIGNAL 

Here we make a remark on how the formulas (0), ( p^ can be applied to the case 
of a non-periodic field. Actually, the expression ( pT]) for the multi-channel characteristic 
function x(A) is quite general and can be applied to the non-periodic case as well. Suppose 
the charge is measured during a finite time T^^\ This can be taken into account by making 
the parameter A of the generating function time-dependent, and by "turning it on" only for 
the time interval of the measurement. If the external field is encoded into a (non-periodic) 
phase factor f(t), then the distribution is given by expressions (|23|) and (|T2|), where 
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A = A(t) = ( J 



:di) 



m 



A, < t < TW; 

0, t<Oort>T(°). 



(D2) 



The operator S(t) becomes unity outside the observation time interval, and this makes the 
determinant ([TT|) well defined. We postpone a general discussion of the non-periodic case 
to elsewhere. Here we only demonstrate that this approach gives the correct answer for the 
equilibrium noise {{n^)) in the absence of external field. 
By the argument of Appendix B, 



-Tr 



92 



9A2 



A=0 



ln(l + fi(S-l)) 







lim 



2\A\ 



dU 



1 



2\A\ 



2 A 



ti-t2 + iO 
2 1^'°' dh 



+ 2\Ar\B 



2m {ti ~ t3 + i0){t3 - t2 + iO) 

dh 1 







-2U|2 



lim 

2'Ki t2-*ti 

dt 



2m (ti-t3 + z0)(t3-t2 + ^0) 

dt3 1 



2m {ti - ts + i0){t3 - t2 + iO) U - t2 + iO 



\A\ 



In 





T(o) 



2m 7,^'<% 27ri (t - ts)^ 27r2 Jo 



rp(0) 



dt 



T(o)-t t 



(D3) 



where r^c is the ultraviolet cutoff set by a characteristic scattering time of the system h^^^. 

It is straightforward to check that the fiuctuations given by the last line of Eq. (P3|) agree 
with the Nyquist equilibrium noise spectrum 5"^^ = ^G{uj)\uj\. Indeed, 



dt I dt'mw))) 



JO 

dco 11 -e*^T 



M|2 2^(0) 

— TT- In 



where \A\ 



\A{E, 



27r UJ^ TT^ 

G'(co' = 0), and h/rgc is the ultraviolet frequency cutoff. 



(D4) 



APPENDIX E: GAUGE TRANSFORMATIONS 



In this Appendix we show how by gauge transformations one can switch between the 
problem with an electric potential and that with a vector potential. Recall that the gauge 
transformation 



changes the vector potential a{x,t) and the electric potential v{x,t) as 



(El) 
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cd(f>(t,x) 

a = a 



e dx 

v = v+'-^-^. (E2) 
e ot 

Assume that we apply a time-dependent magnetic field with the vector potential a{x, t) 
localized around the scatterer. For example, this can be realized as a varying magnetic fiux 
threading a conducting loop with the contact. Then, by the gauge transformation ( [El| ) with 



e 

(j){t,x) = - a{t,x)dx (E3) 

C J —oo 

we can turn to the problem with the zero vector potential and the electric potential 

= (E4) 
In this case, the gauge phase shift across the scatterer is 



e 



+00 



^(t) = <j){t, x) llztZ = - / a{t, x) dx, (E5) 



c 

and it can be viewed as the time-dependent voltage 

applied to the contact. 

Thus we recalled the familiar gauge transformation for the one-particle problem. When 
looking at a many-body problem, one also has to transform the density matrix and Green's 
functions according to the rule (^). Let us consider the occupation number operator of a 
reservoir: 

nH=^5(^-E,)^i4, (E7) 

k 

where the summation is over all reservoir states. In the determinant we have n = (fi) 
averaged over actual distribution. In the time representation 

n(t, t') = E e^^'=(*-*')(^i^.) = E e^^'=(*-*')np(E.) , (E8) 

k k 

where ripiE) is Fermi distribution Under the gauge transformation ( [El|) with the gauge 
phase ( |E3| ) the occupation number operator transforms as 

fiL(t,t') = nL(t,t') 

hnit,t') = nnit,ty^^^'^-^^''^^ (E9) 

for the left and the right reservoirs respectively. 

We assume that the scattering is instant, then the alternating field effect on the scattering 
states is entirely determined by the phase shift v?(t). This means that the ac magnetic field 
can be introduced simply by adding the phase in the scattering amplitudes , while in 
the problem with the electric field we have the energies of occupied states in one reservoir 
to be shifted with respect to those in the other one. By virtue of gauge invariance, these 
two formulations are obviously equivalent, and we make use of it in the discussion of the 
determinant regularization (see Eq.(p3|)). 
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Figure Captions 



Figure 1 : Voltage time dependence (@), (^Sj) , and (|29|) that makes current fluctuations 
minimal. Shown is the simplest solution: periodic sequence of Lorentzian peaks of the area 
h/e each. For the optimal pulse shape the fluctuations do not depend on the pulse width r. 



Figure 2 : Voltage optimal time dependence corresponding to a pair of 27r— kinks of ^pit) = 
f /-oo V{t')dt'. For the pulses of equal sign (a) the counting statistics do not depend on the 
relative position of the pulses ti,2 and on their duration ri 2 (see Example 2, Eq.(|55D). For 
the pulses of opposite sign (b) average current is zero and other statistics show non-trivial 
interference of the kinks (see Example 3, Eqs.(^)), (|40|)). 
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